Astronomy & Astrophysics manuscript no. paper 


©ESO 2011 


February 1,2011 





Alleviating a quenching by solar wind and meridional flows 

Dhrubaditya Mitra 1 , David Moss 2 , Reza Tavakol 3 , and Axel Brandenburg 1 ' 4 



NORDITA, AlbaNova University Center, Roslagstullsbacken 23, SE-10691 Stockholm, Sweden 
School of Mathematics, University of Manchester, Oxford Road, Manchester M13 9PL, UK 

Astronomy Unit, School of Mathematical Sciences, Queen Mary University of London, Mile End Road, London El 4NS, UK 
Department of Astronomy, Stockholm University, SE 10691 Stockholm, Sweden 



February 1,2011, Revision: 1.90 



o 

(N 
C 

(N 



ABSTRACT 



Aims. We study the ability of magnetic helicity expulsion to alleviate catastrophic a-quenching in mean field dynamos in two- 
dimensional spherical wedge domains. 

Methods. Motivated by the physical state of the outer regions of the Sun, we consider a 2 Q mean field models with a dynamical a 
quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow 
("solar wind") and meridional circulation. 

Results. We find that a wind alone can prevent catastrophic quenching, with the field saturating at finite amplitude. In certain param- 
eter ranges, the presence of a large-scale meridional circulation can reinforce this alleviation. However, the saturated field strengths 
are typically below the equipartition field strength. We discuss possible mechanisms that might increase the saturated field. 

Key words. Sun: dynamo - magnet ohydrodynamics (MHD) 
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43 1 . Introduction 

Q • Mean field dynamo models have provided an important 
1 framework for studying the generation of large-scale astro- 

c/3 \ physical magnetic fields and their spatio-temporal dynam- 
. ics. However, these widely used models have been pre- 
sented with a serious c hallenge - namely the so call ed catas- 
, trophic a quenching dGruzinov & Diamond Il994l) . In the 

■ mean field (MF) context this effect, which is a consequence 
I of th e conservation of magnetic helic i ty (|Krause & Radlen . 

(N '■ 119801 IZeldovich. Ruzmaikin & Sokoloffl Il983h . manifests it- 
' self as the decrease of the q -effect with increasing mag- 
nI ' netic Reynolds number R eM dVainshtein & Cattaneol 1992; 
OO ■ ICattaneo & HughesL 1 19961) at finite field strength. In models 
f^*) ' without magnetic helicity fluxes, the quenching of a can become 
<0 severe, with a decreasing as Re^ 1 - truly catastrophic for dy- 
^ — 1 ■ namo action in the Sun, stars and galaxies where the Reynolds 
I numbers are all very large (> 10 9 ). This catastrophic quench- 
• i-H , ing is captured by mean-field models which use dynamical alpha 

■ quench ing, such as that considered by Blackman & Brandenburg 
H \ (120021) . This catastrophic quenching is independent of the de- 

. . . i tails of the dynamo mechanism and is a direct effect of the con- 
servat ion of magnetic helicity, see, e.g.. lBrandenburg & Kapvla 
(120071) who have de monstrated catast r ophic quenching for non- 
local alpha effect or IChatteriee et alj d201 If ) who have demon- 
strated the occurrence of catastrophic quenching in distributed 
dynamos. It has been suggested that the quenching may be al- 
leviated by the expulsion of magnetic helicity through open 
boundaries dBlackman & Field 120001: iKleeorin et alll2000h . At 
least three different physical mechanisms may help in the ex- 
pulsion of small s cale magnetic helicity: (a) large scal e shear 
( Vishniac & Cho, 2001; Subramanian & Brandenburg, 2004; 
Brandenburg & SandinL l2004t iMoss & Sokoloffl. 1201 ll): (b) tur- 
bulent diffusion of magnetic helicity dMitra et all 1201 Oab : (c) 
non-zero mean flow out from a boundary of the domain, e.g. 



a wind. A number of recent studies have demonst rated the possi- 
bility of this alleviation of quenching for solar |Chatteriee et all 
1201 ll 120101: iGuerrero Chatterjee. & B randenburg. 2010) and 
galactic dynamos (e.g. lShukurov et afl.l2006l) . 

In this paper we study the effects of a number of mecha- 
nisms which may facilitate the expulsion of magnetic helicity 
from the dynamo region. Initially we consider t he effects of ad- 
vectio n by a mean flow in a similar manner to Shu kurov et al.l 
d2006l); see also the recen t study in a one dimensional model by 
iBrandenburg et all d2009l) . We envisage that in the Sun the wind 
could be loaded with magnetic helicity th rough coronal mass 
ejections dBlackman & Brandenburgll2003l) . Another potentially 
important mechanism is meridional circulation. The presence of 
such a circulation in the Sun is supported by a number of obser- 
vations which have found evidence for a near-surface poleward 
flow of 10 — 20 ms _1 . Even though the corresponding compen- 
sating equatorward flow has not yet been detected, it is however 
assumed it must exist because of mass conservation. Substantial 
effort has recently gone into the construction of flux transport 
dynamo models which differ from the usual a£l dynamos by 
including an additional advective transport of magnetic fl ux by 
meridional circulation, (see e.g. Dikpati & Gilman, 2009i for a 
recent summary). If magnetic flux is advected by meridional cir- 
culation, it can be expected that such a circulation will also trans- 
port magnetic helicity to the surface layers, which might thus fa- 
cilitate its subsequent expulsion by the wind. We therefore study 
the effects of meridional circulation on the quenching. 

The structure of the paper is as follows. In Section [2] we in- 
troduce our model and its various ingredients. Section[3]contains 
our results, and we give a short summary here. First we con- 
sider our model with an imposed wind but no meridional circu- 
lation. We show that a strong enough wind that penetrates deeply 
enough into the convection zone can indeed alleviate quenching. 
We then make a systematic study of the alleviation of quenching 
as a function of the two parameters specifying the wind, namely 
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the maximum velocity and the depth down to which the wind 
penetrates the convection zone. Next we select a particular set of 
these two parameters such that for large ReM there is no allevi- 
ation of quenching. We then introduce a meridional circulation 
and show that a combination of a wind and circulation is able to 
limit the quenching in cases where the wind alone cannot. We 
further study the effect of the characteristic velocity of merid- 
ional circulation on quenching. Our conclusions are presented in 
Section|4] 



2. The model 

We study two-dimensional (axisymmetric) mean field models 
in a spherical wedge domain, n < r < r%, &\ < 9 < ir/2, 
where r, 9, <fi are spherical polar coordinates. The choice of this 
"wedge" shaped domain is motivated by recent Direct Numerical 
Simulations (DNSs) of forced and connective dynamos in spher- 
ical w edges cut from spheres (Mi tra et al.l 12010b; Kapvla et al.L 
1201 Oh . and the intention to make a similar development of this 
work. 

We consider an a 2 ft mean field model with a "dynamical 
alpha" in the presence of an additional mean flow U. In the sim- 
plest case, where we consider no wind and no meridional circu- 
lation, the mean flow is in the form of a uniform rotational shear 
given by U = [/ shcar = <f>S{r — rrj) sin 9. For the more realistic 
cases we use 



U = U ; 



shear 



u 



wind 



u c 



(i) 



where jy wind and [/ Clrc are respectively the large-scale veloc- 
ity of the wind and circulation. The particular forms we use are 
given in Sections [2. ll and l2~2l below. Thus, we integrate 



d t B = V x (17 xB + S) + r]V 2 B, 



d t a M 



-2rjki 



SB 

B 2 



Re 



M 



V • (a u U), 



(2) 
(3) 



where £ = aB — %J, and a = au + ax is the sum of 
the magnetic and kinetic a-effects respectively. The magnetic 
Reynolds number, Rejvt/3 = f]t/r) and B eq is the equipartition 
field strength. We take rj t — 1, B cq — 1 and fcf = 100 in our 
simulations. Here Eq. (0 is the standard induction equation for 
mean fiel d models and Eq. (f3]> describes the dynamical evolution 
of a; see Blackman & Brandenburg (2002). The last term in the 
right hand side of Eq. (O models the advective flux of magnetic 
helicity. 

We solve Equations © and (0 using the PENCIL Codi£] 
which employs a sixth order centered finite-difference method 
to evaluate the spatial derivatives and a third order Runge-Kutta 
scheme for time evolution. 

Our aim here is to study the effects of the various mecha- 
nisms discussed above in alleviating the catastrophic quenching 
of the magnetic field as ReM increases. 

2. 1. The wind and the "corona" 

In order to include the effects of the solar wind we must include 
an outer region in our model through which the wind flows, by 
extending the outer boundary beyond the convection zone to ra- 
dius r% > r 2 . We shall refer to the region r 2 < r < r% as the 
'corona' . We take the wind to be strong in the corona and to grow 
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weaker as we go into the convection zone. This is represented by 
choosing the following form for [/ Wlnd ? 



U, 



wind 



u. 



wind 



wind 



o, 
o. 



1 + tanh 



''2 



W 



(4) 

(5) 
(6) 



where Uq and w are control parameters which determine the 
strength of the wind speed and its depth of penetration into the 
convection zone respectively. Larger values of w correspond to 
deeper penetration. We let the kinematic a-effect to go to zero 
in the corona by choosing 

-tt/2 



ot\c = tann 

2 

with ocq = 16. 



0.05 



1 — tanh 



r - r 2 



(7) 



2.2. The meridional circulation 

We consider the effects of a meridional circulation, by including 
a velocity U cnc given by 

1 d 



JJCVC = 



-'amp 



9(r) 



US 1IC = 



sin 9 89 



(sin 9 ip) . 



f(r) 

g( r ) 



r dr 
(0-0i)cot< 



(r 
1 
2 
0. 



1 — tanh 



r 2 



U'c 



(8) 

(9) 

(10) 
(11) 
(12) 

(13) 



Here w amp is a parameter controlling the magnitude of circula- 
tion speed and w c i IC determines the effective depth of penetra- 
tion of the circulation into r > r 2 . As a characteristic speed of 
circulation, v C i IC , we take the maximum absolute magnitude of 
the 9 component of U C1TC at r = r%, i.e. at the surface of the 
Sun. Helioseismology shows this velocity to be about 10 to 20 
metres per second in the Sun. A typical velocity field is shown in 
Fig-E and the profile of «k is shown in Fig. [2] In these Figures, 
the parameters are chosen to be a>o — 16, w a — 0.2, Uq = 2, 
r 2 = 1.5, w = 0.3, u amp = 75, r circ = 0.98, w circ = 0.02. 

2.3. Boundary conditions 

For the magnetic field we use perfect conductor boundary con- 
ditions both at the base of the convection zone (at r = r{) and 
at the lateral boundary at the higher latitude (9 — 9\). We as- 
sume the magnetic field to be antisymmetric about the equa- 
tor (9 — tt/2), and at the outer radial boundary of the corona 
(r = r$) we use the normal field condition. In terms of the mag- 
netic vector potential A = (A r , Ag,A$), where B — V x A, 
these conditions become 
f)A 

' " 1 ' (r = n), (14) 



dr 

A r = 0, 



A e = A 4 
dA g 



An 



A r -- 
dA r 



8Ag 

89 
= A fi 



A^ 

d_ 

09 



dA^ 
dr 

= 0x), 



Art 



(r = r 3 ), (15) 



(sin6A <t> )=Q (6 = tt/2). 



(16) 



(17) 
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Fig. 3. Spacetime diagrams for the three components of the mag- 
netic field. These plots are for the set of parameters correspond- 
ing to the velocity field and q;k shown in Figs.[T]and [2] 



Fig. 1. Plot of the velocity field: the arrows show the meridional 
circulation and the wind, and the contours show the angular ve- 
locity. The solar radius is taken to be unity. Although our domain 
extends out to 5 solar radii, for clarity only a part of it is shown 
here. The curve at unit radius denotes the surface of the Sun. 




Fig. 2. The kinetic alpha effect, ok, and wind radial velocity, U r , 
as a function of radial coordinate r for three different latitudes, 
equator (upper curve), mid-latitude (middle curve) and latitude 
of upper boundary (lower curve). Note that the curves for the 
radial velocities differ only in r < 1, where the meridional cir- 
culation is non-zero. 



For «m, on those boundaries where the boundary condition on 
the magnetic field is "perfect-conductor" (i.e. at the bottom of 
the convection zone and at the higher latitude), we choose 



a M = 0. 



(18) 



At the other two boundaries, we recall that since the PDE being 
solved is of first order in space we only need to specify one con- 
dition, which we have already imposed at the lower boundary. To 
calculate the derivative at the outer boundary we therefore just 
extrapolate the solution from inside to outside by a second or- 



der polynomial extrapolation. This is equivalent to using second 
order one sided finite difference at these boundaries. 

As the initial condition for the magnetic field we choose our 
seed magnetic vector potential from a random Gaussian distri- 
bution with no spatial correlation and root-mean-square value 
of the order of times the equipartition field strength. Also, 
initially we take au = q - «k = 0. 



3. Results 

In order to demonstrate that our dynamo is excited, and dis- 
plays both oscillations and equatorward migration, we first use 
the velocity field and the kinetic a profile shown in Fig.[T| with 
ReM = 3 x 10 2 and solve Eqs. (f2j) and (0 simultaneously. The 
resulting space-time diagram for the three components of the 
magnetic field is shown in Fig. [3] This is a typical example of 
the "butterfly" diagrams that are obtained with this model. 

As mentioned above, an important feature of MF dynamos 
in the absence of wind and meridional circulation (i.e. when 
Uo — and v amp = 0), is that they are severely quenched 
as ReM increases. To show this we have plotted in Fig. [4] (a) 

the time-series of the total magnetic energy Em = h (B ) for 
several values of ReM- Here, (...) denotes averaging over the 
domain r\ < r < ri. Clearly the total magnetic energy de- 
creases with ReM- Similar quenching, as a result of the dynam- 
ical evolution of the alpha term, has been s een in many dif- 
feren t mode ls of the solar dynamo (see, e.g.. IChatteriee et ail 
120111 120101: iGuerrero. Chatteriee. & Brandenburg! 120101 for 
some recent examples) , and also in models of galactic dynamos 
(IShukurovetal.Ll2006h . 

To substantiate this further we plot in Fig. the time- 
averaged magnetic energy (Ey[)t as a function of ReM, where 
the time averaging is done over several diffusion times (T) in the 
saturated nonlinear stage (i.e. after the kinematic growth phase 
is over). Time averaging is here indicated by the subscript t after 
the averaging sign. As can be seen in the absence of wind, i.e. 
with Uo = 0, such time-averaged energy falls off approximately 
as Re^ 1 . This gives a quantitative measure of the quenching. 
(The point at ReM = 2 x 10 6 appears anomalous; we believe this 
is because we have not run the code for long enough to achieve 
the final saturated state.) 
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Fig. 4. (a) Time series of Em for Uq — and i; amp = 0, i.e. 
no wind or circulation, for 5 different values of ReM- (b) The 
corresponding plot with f/rj = 1, W = 0.3 and i; amp = 75. 
The time-averaged magnetic energy (i?A/)i is calculated by 
time-averaging these time-series between the two dashed ver- 
tical lines. The other parameters used are otQ = 16, w a = 0.2, 
r 2 = 1.5, and w = 0.3. 



To demonstrate the ability of the wind and circulation to act 
together to alleviate quenching, we have plotted in Fig. [4] (b) 
the time-series of Em for several different values of ReM, in 
the presence of the wind (with Uq = 1, "amp = 75 and depth 
parameter w — 0.3). The dependence of the time-averaged mag- 
netic energy (Em)i on ReM in this case is also plotted in Fig. [5] 
Comparing Fig.|4](b) with Fig.|4](a) and also comparing the two 
lines in Fig.[5]we clearly see that with the parameters chosen the 
wind in conjunction with the circulation is capable of alleviat- 
ing quenching. This is one of our principal results. Note that the 
saturated mean field energy that we observe at large ReM is still 
rather small, only slightly exceeding 10~ 4 of the equipartition 
value. 

Next we attempt to isolate the role of each parameter in our 
model. First we make a detailed systematic study of how quench- 
ing depends on the two parameters Uq and w of our model, for 
a fixed value of ReM = 10 7 and zero circulation, v amp = 0. For 
each pair of parameters we ran our code for up to 50 diffusion 
times. In some cases the time series of Em declines as a function 
of time initially, but at larger times recovers to unquenched val- 
ues, e.g. Uo — 1 in Fig. [6] In some other cases we observe that 
the recovery is merely temporary and at large times Em goes 
to zero. As an example we first show in Fig. [6] the time-series 
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Fig. 5. The time-averaged magnetic energy as a function of ReM 
for Uq = (no wind) and Uq = 1 and w = 0.3 and v amp = 75. 



of Em for various values of Uq, for a fixed w = 0.3. Clearly, 
as the wind velocity increases the transport of magnetic helicity 
out of the domain at first becomes more efficient and we observe 
less quenching. But this alleviation of quenching must have its 
limits because for a large enough wind speed the magnetic field 
itself will be advected out of the domain faster than it is gener- 
ated, thus killing the dynamo (see, e.g., IShukurov et ail [2006 ; 
iBrandenburg et all 119931; iMoss et al.L l2010b . However with pen- 
etration factor w = 0.3 we did not find this effect, even when 
Uq = 100, but with w = 0.5, winds with Uq > 20 kill the dy- 
namo. We deduce that it is necessary to advect large-scale field 
from a substantial proportion of the dynamo region for the dy- 
namo to be killed by advection. 

Then we consider the parameter w which controls the depth 
of penetration of the wind into the convection zone. The depen- 
dence of the time-series of magnetic energy on this parameter is 
shown in Fig. [7] for ReM = 10 7 and Uq — 2. We also note that 
there is a subset of parameters for which the transients are so 
long that it is difficult to decide whether the asymptotic state is 
a quenched dynamo or not, within reasonable integration times. 
In our parameter space, i.e. in the Uq — w plane, the positions 
of the quenched and unquenched runs are shown in Fig.|8j sum- 
marizing the dependence of quenching on these parameters. For 
all the runs we label as unquenched the butterfly diagram is also 
restored at large times. 

3. 1. The effect of circulation 

Next we consider the effect of meridional circulation on the 
quenching. If the wind penetrates inside the convection zone too 
deeply then we expect that circulation will have either no effect, 
or just a marginal effect, because the wind by itself will be ef- 
ficient enough in removing small-scale magnetic helicity from 
deep within the domain. But if the wind does not penetrate so 
deeply, circulation may play an important role in dredging mag- 
netic helicity from deep inside the domain to near the surface 
from where the wind can remove it. To see whether this idea 
can work, we select one point in the phase diagram in Fig [8] 
where we obtain the quenched solution marked by the arrow. 
Then we turn on the meridional circulation. The comparison 
between the time-series of Em with and without circulation is 
shown in Fig. [9] It can be seen that the final magnetic energy 
reached does not depend on the amplitude of circulation if the 
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Fig. 6. Time series of Em for 5 different values of Uo, namely, 
U Q = 0.1(*), 0.4(0), 1(A)2(D) and 10(x) with all other param- 
eters held fixed, in particular ReM = 10 7 , w = 0.3. 
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Fig. 7. Time series of Em for 5 different values of w, [w — 
0.01(*),0.2(A), 0.3(D), and 0.4(x).] with all other parameters 
held fixed, i.e., ReM = 10 7 , Uq = 2 with no circulation. 
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Fig. 8. The incidence of quenched and unquenched solutions in 
the Uq—w plane. The positions of quenched models are denoted 
by the symbol *, the symbol identifies unquenched models. 
The arrow is explained in Sect. I3.ll 



Fig. 9. Time series of Em as a function of time for Uq = 2 and 
w = 0.25 with w amp = 0(+), 75(*), 100(0), and200(A). 

amplitude of circulation is greater than a critical value. Note that 
this alleviation of quenching by the circulation only works for 
those points in the Uq — w parameter space which lie close to the 
boundary between the quenched and non-quenched states in the 
phase diagram. For points with very small w, i.e. in cases where 
the wind penetrates very little into the convection zone, even a 
very strong circulation cannot remove the quenching. 

Another possible mechanism that can transport magnetic he- 
licity from the bulk of the convection zone to its surface is 
the diffusion of magnetic helicity. This can be described by 
adding the term K t V 2 aM to the right hand side of Eq. (01, where 
K t is an effective turbulent diffusivity of the magnet ic helicity . 
Nume rical simulations have estimated K t ~ 0.3?7t (iMitra et aU 
2010a). We have checked that such a diffusive flux of mag- 
netic helicity can alleviate quenching at least as effectively as 
the meridional circulation, in the presence of the wind. 

Finally we note that the alleviation of quenching as described 
here is independent of some details of the underlying dynamo 
model. In particular it does not depend on whether we have an 
a 2 dynamo or an a 2 il dynamo. To check this assertion explic- 
itly we also solved the same problem but with U shBaI = in 
Eq. (JTJ. The results are shown in Fig. [TOl where we compare the 
alleviation of quenching for the a 2 dynamo (top panel) against 
the corresponding a 2 Vl dynamo (bottom panel). 



4. Conclusions 

We have introduced two observationally motivated effects that 
may help reduce the catastrophic quenching found in mean 
field dynamo models. An outward flow from the dynamo region 
("wind") is found to be effective in allowing the quenching to 
saturate at finite values of the field strength. The wind alone is, 
however, only effective when it penetrates quite deeply into the 
convection zone. These effects are modified to some extent by 
the presence of a meridional circulation which has the ability to 
transport small scale helicity from deep in the convection zone 
to near the surface, from where the wind can more effectively re- 
move it. However, the effects of circulation in our model are not 
dramatic. It is also true that the saturation fields in our model are 
rather small compared to the equi partition field strength. This 
was also observed in t he mo del of IShukurov et al.l (120061) : see 
also lMoss & Sokoloffl (1201 ll) . One possibility, that we have not 
explored, is that the neglected inhomogeneity of the solar con- 
vection zone may be important. 
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Fig. 10. The behaviour of the time-averaged magnetic energy as 
a function of magnetic Reynolds number ReM, which shows the 
alleviation of quenching, with wind speed Uq = 2 and depth pa- 
rameter w = 0.8 (upper panel). Also shown is the corresponding 
plot in the absence of a wind, which clearly shows a catastrophic 
quenching. The lower panel is for a a 2 il dynamo. 



It is interesting to try to estimate the various parameters of 
our model in physical units and to compare them with the so- 
lar values. We have taken the solar radius as the unit of length 
(7 x 10 10 cm). The a effect can be taken to be a measure 
of the small scale velocity in the Sun, a ~ (1/3) where 
u = U — U. The iBaker & Temesvarvl (Q%6) tables give es- 
timates for small-scale velocities in the convection zone of the 
Sun of between 4 x 10 3 -2 x 10 5 cm s _1 , in regions where con- 
vection is efficient. As we have considered the convection zone 
to be homogeneous we consider 10 4 cm s^ 1 to be a reasonable 
estimate. Then, as a = 16 in our units the unit of velocity is 
- 10 4 /(3 x 16) cms -1 — 2 x 10 3 cms -1 , and the unit of time, 
obtained from length and velocity units given above, is ~ 10 s 
s m 10 yrs. Thus our characteristic cycle period, T m 1, cor- 
responds to approximately 10 years. Then the maximum wind 
speed we have used (Uq = 10) would correspond to 2 x 10 3 
cm s _1 . The speed of the meridional circulation at the surface in 
our units is v SUI f — 0.47 for v amp = 75. Translated to physical 
units this becomes v SUI { « 1 m s _1 , which is of the same order 
of magnitude as the solar meridional velocity. If in the estimates 
above we use the maximum and minimum values of the small- 
scale velocity as given by the Baker & Temesvary tables, instead 
of the mean, the maximum surface speed of meridional circula- 
tion will be between 0.4 m s -1 and 20m s _1 . The speed of the 
solar wind that we have used is significantly smaller than that of 
the actual solar wind, but on the other hand the real solar wind is 



a highly fluctuating turbulent flow, whereas we have considered 
a constant outflow. 

To summarise, we have presented a very simplified model, in 
order to explore some basic ideas relevant to the solar dynamo. 
We cannot claim to have "solved" the quenching problem, but 
feel we have identified, and to some extent quantified, mecha- 
nisms of potential interest. We appreciate that there are a number 
of desirable improvements, even in this MF formulation. These 
include using a more realistic solar-like rotation law, investiga- 
tion and comp arison of the effe cts of other fluxes of magnetic 
helicity (e.g. Zhan g^staLj, | 2006l) . the diffusive magnetic helic- 
ity flux dMitra et all 1201 Oal) . the inclusion of compressibility in 
some form, but most importantly perhaps, using a more realistic 
model for the solar wind allowing for magnetic helicity load- 
ing via coronal mass ejections. Notwithstanding these possible 
shortcomings, we do feel that our results provide motivation for 
further investigations in the co ntext of solar and stellar dynamos . 
Investigations using DNS (e.g. lWarnecke & Brandenburg, 2010) 
appear likely to be especially interesting, and we hope to pursue 
this approach. 
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